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ABSTRACT 

(Distribution  Limitation  Statement  No.  3) 


Gas  gun  impact  experiments  were  performed  on  five  thicknesses  of  6061-T6  aluminum 
targets  to  determine  the  incipient  spall  threshold.  Data  obtained  were  compared 
with  previous  data  for  the  same  material  obtained  by  General  Motors  and  other 
laboratories.  Results  indicated  that  batch-to-batch  variations  in  properties  of 
a  commercially  "standard"  material  may  cause  up  to  15  percent  differences  in  spall 
thresholds  for  short  duration  (<0.1  u  second)  shock  loads.  Four  mathematical 
models  of  dynamic  fracture  were  evaluated  with  the  spall  lata  obtained.  The 
Hole-Growth  model  developed  by  Stanford  Research  Institu.e  under  contract  to  the 
Air  Force  Weapons  Laboratory  proved  superior  to  the  other  models  considered. 

This  superiority  was  manifested  in  the  ease  of  interpretation  of  the  calculated 
damage  levels  and  the  physical  significance  of  the  calculations. 
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SECTION  I 
INTRODUCTION 

Studies  of  dynamic  failure  in  homogeneous  materials  have  been  of  particular 
interest  to  the  USAF  and  DOD  because  of  the  applicability  of  such  studies  to 
weapons  systems  subjected  to  a  nuclear  detonation.  Absorption  of  a  given  dose 
of  X  rays  by  a  weapons  system  causes  a  high-amplitude,  short-duration,  compres¬ 
sive  shock  pulse  to  propagate  through  weapon  heat  shields  and  metal  substructures. 
This  pulse  is  reflected  from  a  free  surface  or  interface  as  a  relief  wave  propa¬ 
gating  in  the  opposite  direction.  The  interaction  of  relief  waves  causes  tension 
in  the  material.  The  behavior  of  the  material  as  a  result  of  this  tension  is 
important  because  catastrophic  failure  or  significant  degradation  of  structural 
properties  may  result. 

Many  studies  have  been  made  of  time-dependent  fracture  in  homogeneous  materials. 
One  material  which  has  undergone  extensive  investigation  is  6061-T6  aluminum 
because  its  behavior  under  dynamic  loading  conditions  appears  to  be  relatively 
strain-rate  insensitive  (figure  1).  The  Air  Force  Weapons  Laboratory  (AFWL)  under¬ 
took  a  study  of  6061-T6  aluminum  to  determine  its  spallation  threshold.  This 
report  provides  a  description  of  the  experimental  technique,  a  summary  of  the 
results  obtained,  a  comparison  of  those  results  with  data  obtained  by  other  labora¬ 
tories,  and  an  evaluation  of  a  few  of  the  mathematical  models  of  dynamic  fracture 
currently  available. 
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Figure  1.  Strain  Rate  Sensitivity  of  6D61-T6  Aluminum  at  Room  Temperatur 
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SECTION  II 

EXPERIMENTAL  TECHNIQUE 

As  mentioned  In  section  I,  absorption  of  X  rays  by  a  materia]  gives  rise  to 
a  compressive  pulse  of  high  amplitude  and  short  duration  which  reflects  from  free 
surfaces  or  material  interfaces  to  produce  tension.  In  the  homogeneous  materials 
such  as  metals,  this  tension  reaulta  in  the  nuclcatlon  of  microcracks  or  voids 
and  their  subsequent  growth.  When  the  tensile  stress  reaches  still  higher  ampli¬ 
tudes  or  longer  durations,  these  cracks  begin  to  coaleece  and  total  separation 
of  the  material  can  occur. 

Slmualtlon  techniques  exist  through  which  comprasslve  shock  wavss  similar  to 
those  caused  by  X-ray  deposition  can  be  Introduced  into  s  materiel.  One  such 
technique  Is  the  hlgh-veloclty  Impact  on  a  target  material  by  s  llyar  plsle. 

This  was  the  method  chosen  for  the  etudy  being  reported.  A  schematic  drawing  of 
the  experimental  apparatus  Is  shown  in  figure  2.  The  flyer  piste  Is  accelerated 
by  a  sabot  to  its  Impact  velocity  which  le  measured  by  three  velocity  pins.  After 
the  flyer  plate  Impacts  the  target,  the  sabot  Is  stopped,  both  target  and  flyer 
are  recovered  from  a  catching  mechanism. 

Figure  3  shows  a  cross-sectional  view  of  both  the  target  and  flyer  plates. 

All  plates  were  mechlned  from  1/4  Inch  thick  plate  stock.  The  ratio  of  target 
thickness  to  flyer  thickness  was  kept  at  Lite  constant  value  of  7:1.  There  were 
two  reasons  for  the  choice  of  this  ratio.  First,  making  the  target  twice  the 
thickness  of  the  flyer,  end  of  the  same  material,  forces  the  duration  of  the  shock 
pulse  to  be  approximately  equal  to  the  transit  time  of  the  shoe1-  across  the  tar¬ 
get  and  thereby  minimizes  attenuation  of  the  puloe  duo  to  the  trail  lug  rarefaction 
wave.  The  choice  of  this  ratio  also  causes  spallation  to  ocou:  at  approximately 
the  midplane  of  the  target.  The  first  Intel  action  of  two  rarefaction  waves  and, 


3 


Figure  2.  Schematic  Drawing  of  Experimental  Apparatus 
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Figure  3.  Target  Plate  Cross  Section 
and  Flyer  Plate  Cross  Section 
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thus,  the  first  tensile  stress,  occurs  at  this  plane.  The  simple  X-t  curve  shown 
in  figure  4  demonstrates  both  of  these  points  and  also  shows  that  the  midplane 
of  the  target  is  under  tensile  stress  for  the  longest  period  of  time.  herefore, 
the  microcracks  that  form  in  the  target  materia?  are  concentrated  around  this 
plane. 

Previous  experimenters  had  demonstrated  the  time-dependency  of  spall  in  homo¬ 
geneous  materials  (Ref  1-6).  The  objective  of  this  experiment  was  to  determine 
the  threshold  for  incipient  spall  in  6061-T6  aluminum.  For  purposes  of  the 
experiment,  incipient  spall  is  defined  as  the  amount  of  damage  that  occurs  such 
that  the  spall  plane  is  50  percent  covered  by  voids  or  cracks.  The  spall  plane 
is  the  plane  which  is  under  tension  for  the  longest  period  of  time. 

In  impact  experiments  such  as  these,  the  duration  of  the  pulse  varies  with 
the  thickness  of  the  flyer,  while  the  amplitude  of  the  pulse  is  proportional  to 
the  impact  velocity.  For  the  original  experiment,  five  different  flyer  thick¬ 
nesses  were  chosen  such  that  the  pulse  durations  would  vary  from  less  than  0.1 
usee  to  greater  than  0.5  usee.  There  would  be  five  flyers  of  each  thickness 
fired  at  a  range  of  velocities  both  above  and  below  the  critical  velocity  for 
incipient  spall  for  that  thickness.  The  end  result  was  to  be  five  data  points 
on  the  incipient  spall  threshold  for  6061-T6  aluminum. 

To  determine  the  degree  of  damage,  the  following  technique  was  employed: 

Each  target  was  recovered  from  the  catcher  mechanism  at  the  completion  of  the 
shot.  Abnormalities  such  as  off-center  impact  or  indications  of  multiple  impacts 
were  noted  at  this  time.  The  samples  were  sectioned,  mounted,  and  polished. 

Each  sample  was  then  examined  for  daraag.  magnifications  of  both  50X  and  100X, 
and  photomicrographs  were  taken.  The  results  of  the  program  are  presented  in 
section  III. 
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Figure  4.  Simple  x-t  Diagram  for  6061-T6  Aluminum 
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SECTION  III 
EXPERIMENTAL  RESULTS 

Each  of  the  recovered  target  specimens  was  examined  for  damage  and  each  was 
classified  according  to  figure  5  on  the  following  page.  The  results  of  the 
first  25  gas-gun  shots  are  tabulated  and  illustrated  in  figure  6. 

Figure  6  shows  the  time-dependence  of  the  spall  phenomenon.  Those  specimens 
which  undergo  tensile  stress  through  a  longer  time  duration  begin  to  fracture  at 
somewhat  lower  stress  levels  than  those  which  are  under  tension  for  a  shorter 
period  of  time. 

This  first  series  of  gas-gun  shots  indicated  some  problems  with  the  experi¬ 
mental  configuration.  Note  in  table  I  that  with  the  thinnest  targets  no  spall 
was  detected  at  all.  However,  with  the  32-mil  targets,  some  damage  was  present 
with  even  the  lowest  velocity  flyer.  The  steep  rise  in  spall  threshold  indicated 
in  figure  6  for  the  thinnest  flyers  (shortest  pulse  durations)  was  not  expected 
and  did  not  compare  at  all  to  previous  data  obtained  for  the  6ame  material. 

Most  notably.  General  Motors,  in  work  performed  under  the  DASA-sponsored  PREDIX 
program,  had  reported  damage  in  20-mil  targets  from  10-mil  flyers  with  a  velocity 
of  approximately  860  ft/sec  (Ref  1). 

Since  the  experimental  technique  used  by  both  AFWL  and  GM  was  practically 
identical,  a  reasonable  explanation  was  sought  for  thi6  discrepancy.  It  is 
known  that  the  gas-gun  impact  technique  provides  a  very  flat  impact.  Planarity 
did  not  appear  to  be  the  problem.  Timing  devices  had  recently  been  calibrated, 
and  measured  velocities  were  accurate  to  ‘2  percent. 

It  was  then  noted  that  all  of  the  specimens  had  some  surface  curvature  and 
that  this  effect  seemed  to  be  more  pronounced  in  the  thinner  targets.  It  was 
possible  that,  at  the  high  velocities  being  used,  the  thin  flyers  were  bowing 
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Figure  5.  Levels  of  Damage  in  Spalled  Target  Specimens 
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Table  I 

RESULTS  OF  SERIES  I  IMPACTS 


Flyer  Thickness 
(mils) 

Target  Thickness 
(mils) 

Velocity 

(ft/sec) 

Results 

10 

20 

1373 

No  spall 

1383 

No  spall 

1455 

No  spall 

1571 

No  spall 

1619 

No  spall 

16 

32 

850 

Not  fired 

983 

Incipient 

1028 

Above  incipient 

1064 

Above  incipient 

1150 

Above  incipient 

20 

40 

864 

Not  recovered 

913 

Incipient 

938 

Above  incipient 

985 

Above  incipient 

1040 

Above  incipient 

AO 

80 

575 

No  spall 

581 

No  spall 

661 

Incipient 

744 

Above  incipient 

790 

Above  incipient 

80 

160 

450 

Not  fired 

490 

No  spall 

534 

Incipient 

560 

Above  Incipient 

650 

Not  fired 

11 
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severely  and  the  actual  impact  velocities  in  the  center  of  the  configuration  were 
much  less  than  chose  being  measured.  This  would  be  the  case  if  the  flyer  or 
target  had  to  be  deformed  before  impact  occurred  as  this  deformation  would  decel¬ 
erate  the  flyer  plate.  Also,  the  thicker  flyers  with  the  same  diameter  would 
have  greater  stiffness,  would  bow  less,  and  would  provide  more  accurate  results. 
Finally,  the  thin  targets  exhibited  some  spallation  near  the  edges  where  the 
measured  velocities  were  probably  accurate,  but  no  spallation  in  the  center  of 
the  targets.  For  this  experiment,  only  observations  near  the  center  of  the  tar¬ 
gets  were  of  interest.  The  results  are  to  be  used  to  validate  models  describing 
dynamic  failure  in  homogeneous  materials.  These  models  would  be  incorporated 
into  one  dimensional,  hydrodynamic  computer  codes  such  as  the  PUFF  computer  code. 

To  observe  failure  which  was  primarily  due  to  one-dimensional  stress  wave  propa¬ 
gation,  only  the  center  region  of  the  targets  was  of  interest.  Near  the  edges, 
two-dimensional  effects  arise  due  to  shear  waves  reflected  from  the  edges. 

A  second  series  of  impact  experiments  was  planned  to  examine  the  effects  of 
bowing  in  the  flyer  plates.  The  results  obtained  for  the  40-mil  flyer  impacts 
in  series  I  compared  favorably  with  incipient  spall  data  for  similar  pulse  dura¬ 
tions  obtained  by  other  laboratories.  Thus  the  diameter-to-thickness  ratio  for 
these  flyers  was  adopted  for  the  thinner  flyers.  At  40  mils,  the  ratio  was 

fl  =  i^OO  .  37.5  (1) 

xf  0.040 

Using  a  ratio  of  40:1  resulted  in  the  diameters  shown  in  table  II. 

The  experimental  matrix  for  series  II  included  five  shots  for  each  flyer 
thickness,  a  total  of  fifteen  shots.  The  velocities  chosen  were  again  varied  so 
as  to  bracket  the  critical  velocity  for  incipient  spall  for  each  flyer  thickness. 
The  results  of  this  series  are  summarized  in  table  III  and  illustrated  in  fig¬ 
ure  7.  Also  shown  for  comparison  in  figure  7  are  the  General  Motors  data  (Ref.  1) . 
The  agreement  between  the  two  sets  of  data  is  better  although  the  AFWL  curve 
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Table  II 

TARGET  AND  FLYER  PLATE  DIAMETERS  FOR  SERIES  II  IMPACTS 


Flyer 

Thickness,  xf 
(mils) 

Target 

Thickness,  xfc 
(mils) 

Flyer 

Diameter,  df 
(in.) 

Target 

Diameter,  dt 
(  in . ) 

10 

20 

0.40 

0.20 

16 

32 

0.60 

0.30 

20 

40 

0.80 

0.50 

Table  III 

RESULTS  OF  SERIES  II  IMPACTS 

Flyer  Thickness 

Target  Thickness 

Velocity 

(mils) 

(mils) 

(ft/sec) 

Results 

10 

20 

663 

No  spall 

726 

No  spall 

836 

No  spall 

970 

Below  incipient 

16 

32 

1214 

Above  incipient 

650 

No  spall 

726 

No  spall 

843 

Below  incipient 

958 

Incipient 

1188 

Above  incipient 

20 

40 

690 

No  spall 

698 

No  spall 

853 

Incipient 

956 

Above  incipient 

1025 

Not  fired 
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still  exhibits  a  sharper  rise  (greater  time  dependence)  for  the  thinnest 
flyers . 

One  final  series  of  impact  shots  was  planned  to  attempt  to  resolve  the  25 
percent  difference  in  spall  thresholds  for  short  pulse  durations.  General 
Motors  laboratory  provided  AFWL  with  some  aluminum  cut  from  the  same  plate  which 
GM  used  for  their  spall  experiments.  The  experimental  configuration  was  identi¬ 
cal  to  the  AFWL  shots,  the  only  variable  being  the  material.  In  this  way  it  was 
hoped  to  determine  .  the  spall  threshold  varied  for  the  same  material  from 
batch  to  batch. 

The  experimental  matrix  chosen  Included  four  flyer  thicknesses.  For  each 
thickness,  two  shots  would  be  fired  at  the  AFWL  critical  velocity  and  two  shots 
at  the  GM  critical  velocity.  The  results  of  this  series  are  summarized  in 
table  IV. 

At  best,  these  results  are  Inconclusive.  For  all  targets  impacted  at  the 
AFWL  critical  velocities,  resulting  damage  ranged  from  incipient  to  severe  spall. 
At  GM's  critical  velocities,  the  20-mil  targets  showed  some  damage  below  the 
incipient  level;  while  results  for  the  32-mll  and  40-mil  targets  were  ambiguous. 
However,  closer  examination  of  these  shots  does  point  out  two  interesting  results. 

In  figure  8  the  incipient  spall  thresholds  determined  by  AFWL  and  by  GM  are 
plotted.  The  results  of  the  series  111  impacts  indicate  a  spall  threshold 
occurring  in  the  band  created  by  previous  data.  In  figures  9  and  10  damage 
resulting  from  Impacts  at  nearly  equal  velocities  in  the  AFWL  aluminum  and  the 
GM  aluminum  is  shown.  It  is  obvious  that  in  this  commercially  "standard" 
material,  the  GM  aluminum  was  more  severely  damaged  than  the  AFWL  aluminum  im¬ 
pacted  at  the  same  velocity. 

Thus,  one  might  conclude  that  difficulties  encountered  in  attempts  to  dupli¬ 
cate  spall  data  from  laboratory  to  laboratory  arc  complex  enough  without  the 
added  confusion  contributed  by  variations  in  batch-to-batch  qualities  of 
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"standard"  material.  In  addition,  Che  spall  threshold  apparently  is  best 
defined  by  a  band  in  which  incipient  spall  will  occur  rather  than  a  single  curve. 


Table  IV 

RESULTS  OF  SERIES  III  IMPACTS 


Flyer  Thickness 
(mils) 

Target  Thickness 
(mils) 

Velocity 

(ft/sec) 

Results 

10 

20 

967 

Below  incipient 

982 

Below  incipient 

1089 

Above  incipient 

1108 

Above  incipient 

16 

32 

841 

No  spall 

849 

Incipient 

938 

Incipient 

950 

Above  incipient 

2C 

40 

782 

Incipient 

783 

No  spall 

857 

Incipient 

864 

Above  incipient 

40 

80 

627 

Below  incipient 

639 

Incipient 

676 

Above  incipient 

1200 
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Figure  8.  Series  III  Data  Compared  with  AFWL  and  GM  Spall  Thresholds 
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SECTION  IV 
DATA  ANALYSIS 

The  purpose  of  these  experiments  was  to  determine  the  time-dependent  spall 
threshold  for  6061-T6  aluminum.  These  data  were  to  be  used  to  validate  vai 
dynamic  fracture  models  proposed  for  inclusion  in  the  PUFF  computer  code, 
is  a  one-dimensional,  Lagrangian,  finite-difference  code  used  to  calculate  the 
response  of  layered  homogeneous  materials  to  high-amplitude,  short-duration  shock 
loading.  At  present,  the  spall  criterion  used  in  the  PUFF  code  is  cricical  ten¬ 
sile  stress  (Ref  7).  This  criterion  is  not  valid  for  metals  which,  in  general, 
exhibit  a  time-dependency  in  their  spall  behavior.  Thus,  a  dynamic  spall  criterion 
would  be  a  valuable  addition  to  the  PUFF  code. 

In  this  section,  some  of  the  dynamic  fracture  models  considered  by  AFWL  for 
inclusion  in  the  PUFF  code  will  be  discussed.  These  models  will  be  compared  with 
AFWL  experimental  data  as  well  as  data  obtained  at  other  laboratories.  Finally, 
the  results  of  including  the  models  in  the  P-PUFF  code  (plate  slap  version  of 
PUFF)  will  be  presented. 

1 .  FRACTURE  MODELS 

A  number  of  dynamic  fracture  models,  both  empirical  and  phenomenological, 
have  been  developed  In  recent  years.  Those  considered  by  AFWL  for  inclusion  in 
the  PUFF  code  are  summarized  ir.  table  V.  Other  available  models  (Ref.  8-11)  were 
not  Included  in  this  study  because  a  similar  model  was  studied,  sufficient  data 
was  not  available  for  proper  analysis,  or  the  model  was  not  easily  adaptable  to 
the  PUFF  code. 

The  cumulative  damage  model  was  developed  by  Tuler  and  Butcher  at  Sandia 
Laboratories.  The  constants  are  evaluated  with  experimental  data  for  a  given 
material  (Ref.  12).  The  term  cQ  has  been  compared  to  both  the  static  fracture 
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stress  and  the  dynamic  yield  strength.  The  exponent  X  can  be  assigned  a  value 
of  1.0  to  correspond  to  an  impulse  criterion  for  failure,  a  value  of  2.0  to 
correspond  to  an  energy  criterion,  or  some  other  value  based  on  a  best  fit  of 
^experimental  data.  Generally,  a  value  is  assumed  for  either  aQ  or  X  and  a  least- 
squares  analysis  is  used  to  solve  for  the  other  two  constants.  The  cumulative 
damage  model  is  independent  of  pulse  shape  since  the  integration  can  be  solved 
for  an  arbitrary  shape. 

» 

Table  V 

SUMMARY  OF  DYNAMIC  FRACTURE  MODELS 


Model _ Equation 


rc 

1. 

Cumulative  Damage 

Jo  (a  -  oo)  dt  “  K 

2. 

Rate  Process 

o  =  A  +  B  log  t 

3. 

Series  Expansion 

0  ”  °o  +  t  +  t‘ 1  +  '  *  * 

4. 

Hole  Growth 

R(t)  *  Rq  exp  [A(o)t] 

The  rate  process  model  of  Cohen  and  Berkowitz  (Ref.  13)  requires  the  evalua¬ 
tion  of  only  two  constants,  which  can  be  done  simply  by  a  least-squares  analysis. 
The  difficulty  with  this  model  is  the  definition  of  the  pulse  duration  t.  For 
a  perfectly  rectangular  pulse,  the  value  is  obvious.  For  an  arbitrary  pulse 
shape,  the  pulse  duration  is  determined  through  consideration  of  impulse  and  a 
pulse  shape  constant  k  defined  in  reference  13. 

The  series  expansion  model  was  developed  by  the  author  for  its  mathematicl 
simplicity  and  the  accuracy  with  which  it  reproduces  experimental  data.  The 
model  has  no  physical  significance,  but  can  use  as  many  terms  as  necessary  to 
obtain  desired  accuracy. 

The  Hole-Growth  model  was  developed  by  Stanford  Research  Institute  (SRI)  unde 
contract  to  AFWL.  The  model  uses  experimental  impact  data  to  determine  the 
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nucleation  race  and  growth  rate  of  voids  or  cracks  as  functions  of  stress  ampli¬ 
tude,  stress  duration,  and  instantaneous  void  size.  Data  for  1145  aluminum 
indicate  that  the  nucleation  rate  is  constant  and  that  the  growth  rate  can  be 
described  by  the  viscous  flow  relation 

R(t)  -  RQ  exp  (A(o)t)  (2) 

Where 

R(t)  «  void  radius  at  time  t 
Ro  =  initial  void  radius 
A(o)  “  growth  parameter 

Both  of  these  functions  are  illustrated  in  figure  11.  Stanford  Research  Insti¬ 
tute  has  written  two  subroutines  for  the  PUFF  code  which  compute  the  initiation 
and  growth  of  voids  in  the  material  and  the  reduced  tensile  stress  due  to  recom¬ 
pression  waves  propagating  outward  from  the  newly  formed  free  surfaces  (Ref.  4). 
Application  of  this  model  to  the  experimental  data  is  limited  because  of  the  lack 
of  data  on  the  fracture  parameters  for  6061-T6  aluminum.  However,  some  PUFF  runs 
vlere  made,  and  these  are  discussed  in  the  final  part  of  this  section. 

2..  COMPARISON  OF  OBSERVED  AND  CALCULATED  SPALL  THRESHOLDS 

To  compare  the  AFWL  incipient  spall  threshold  data  with  the  predictions  of 
the  models  and  with  data  generated  at  other  laboratories,  the  P-PUFF  computer 
code  was  used  to  convert  the  velocity-flyer  thickness  data  to  stress-time  data. 

The  version  of  P-PUFF  used  is  described  in  reference  14,  and  is  basically  an 
elastic-plastic  formulation  using  the  Von  Mises  stress-dependent  yield  criterion. 
The  equation  of  state  data  used  in  the  code  calculations  are  provided  in  table  VI. 
Also  shown  for  comparison  in  table  VI  is  the  equation  of  state  data  used  by  other 
laboratories  for  6061-T6  aluminum  and  data  for  mylar  which  are  discussed  later 
in  this  section. 

The  PUFF  code  was  run  for  each  flyer  thickness  and  the  corresponding  critical 
velocity  for  incipient  spall.  The  existing  fracture  routine  was  suppressed  for 
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these  runs.  A  plot  of  the  stress-time  history  at  the  midplane  of  the  target  was 
obtained,  a  sample  of  which  is  shown  in  figure  12. 

Using  the  printed  code  output  and  the  stress-time  histories  for  the  AFWL 
equation  of  state,  a  peak  tensile  stress  and  pulse  duration  were  determined  for 
each  flyer  thickness,  providing  five  de*-a  points  on  the  stress-pulse  duration 
spa.ll  threshold  (fig.  13).  The  pulse  duration  was  taken  as  the  width  of  the 
pulse  at  4.25  kilobars  according  to  the  convention  established  by  the  PREDIX 
group  (Ref.  19).  Other  laboratories  have  used  the  width  at  peak  stress,  at  zero 
stress,  and  at  half  maximum.  Some  convention  is  necessary  since  these  histories 
are  quasi-retangular  and  the  pulse  width  may  vary  considerably  from  peak  to  zero 
stress.  The  value  of  4.25  kilobars  (the  static  yield  limit)  emerged  from  the 
PREDIX  studies  as  the  value  of  oo  in  the  Tuler-Butcher  spall  model,  and  pulse 
durations  were  taken  at  this  stress  level.  The  indicated  error  bars  of  approxi¬ 
mately  10  percent  in  figure  13  are  based  upon  increments  in  velocities  and  visual 
approximations  of  the  level  of  damage. 

These  data  were  then  used  to  determine  the  constants  for  the  first  three 
models  discussed  in  the  previous  section  which  lend  themselves  to  this  type  of 
analysis.  With  the  Tuler-Butcher  model,  three  calculations  were  performed.  The 
value  of  A  was  fixed  at  1.0  and  a  least-squares  fit  was  used  to  determine  o0  and 
K.  This  was  repeated  with  A  =  2.0.  Finally  oQ  was  set  at  4.25  kilobars  and 
values  were  found  for  A  and  K.  For  the  Rate  Process  model,  a  least-squares  fit 
was  used  to  find  the  two  constants  required.  The  Series  Expansion  model  was 
limited  to  two  time-dependent  terms , 

o  =  o0+£  +  f2  (3) 

and  a  least-squares  fit  was  used  to  determine  values  for  0o ,  A,  and  B.  The  per¬ 
tinent  equations  used  in  the  calculations  are  presented  in  appendix  I.  The 
results  of  these  calculations  are  shown  in  table  VII  and  figures  14  through  16. 
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Figure  12.  Stress-Time  History  at  Spall  Plane  (AFWL) 


Figure  13.  Incipient  Spall  Threshold  for  6061-T6  Aluminum 
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Table  VII 


CONSTANTS 

IN  SPALL  MODELS 

Model* 

Assumption 

Variance** 

Tuler-Butcher 

a.  (a  -  9 .000) t  =  1.134 

A  =  1.0 

0.8098 

b.  (o  -  2.995)2t  =  5.504 

A  =  2.0 

0.2227 

c.  (o  -  4 . 25)  1  *  8 1 1  =  4.521 

a  -  4.25 
o 

0.2496 

Rate  Process 

a.  a  =  6.433  -  13.753  logt 

0.4575 

Series  Expansion 

,  .  2.114  0.071 

a.  a  *  6.630  +  — - —  -  - £ — 

t 

0.0681 

*Unit  of  stress  is  kilobars  and  unic  of  time  is  microseconds.  Units  of  constants 
are  consistent  with  these. 

**The  variance  is  a  measure  of  the  goodness-of-f it .  It  is  the  sum  of  the  squares 
of  the  deviations  averaged  over  the  number  of  data  points. 

None  of  the  models  provides  a  bad  fit  of  the  experimental  data.  As  Ferdman 
and  Jajosky  showed  (Kef.  3)  the  Tuler-Butcher  model  provides  a  better  fit  with 
A  =  2  than  with  A  «  1.  Also,  for  fixed  o0 ,  the  best  fit  value  for  A  is  closer  to 
two  than  to  one.  The  Tuler-Butcher  model  provides  a  somewhat  better  fit  than  the 
Rate  Process  model.  The  Series  Expansion  model  obviously  provides  the  best  fit, 
although  the  negative  value  of  the  constant  B  leads  to  an  unrealistic  behavior  of 
the  spall  threshold  at  very  short  pulse  duration. 

The  next  step  in  the  evaluation  of  these  models  was  to  include  more  of  the 
spall  dat;"1  "’vailable  for  6061-T6  aluminum.  The  data  used  are  summarized  in 
table  VIII. 

Some  discussion  of  these  data  is  appropriate  before  applying  any  of  the  models. 
The  data  from  GM  have  been  mentioned  previously.  It  was  obtained  under  the  PREDIX 
program  using  the  same  experimental  technique  as  the  AFWL  used.  The  Sandia 
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*Mylar  flyers. 


DATA 


Peak 

Stress 

(kb) 


Pulse 

Duration 


Pulse 

Shape 


21.70 

19.08 

16.91 

12.29 
10.14 

17.26 

14.35 
12.84 
11.43 
10.54 

12.58 

10.60 

8.82 

19.30 
9.20 

17.93 

12.11 

8.92 

30.80 

26.19 

21.75 

18.52 

16.48 

10.37 

13.36 


0.0824 

0.1289 

0.1598 

0.3133 

0.6126 

0.0774 

0.1461 

0.4056 

0.7470 

1.0890 

0.5395 

0.8216 

1.7241 

0.0894 

0.1892 

0.2956 

0.6957 

1.3514 

0.0566 

0.1254 

0.1370 

0.0826 

0.1865 

0.6849 

0.2926 


Rectangular 

Rectangular 

Rectangular 

Rectangular 

Rectangular 

Rectangular 

Rectangular 

Rectangular 

Rectangular 

Rectangular 

Rectangular 

Rectangular 

Rectangular 

Triangular 

Triangular 

Rectangular 

Rectangular 

Rectangular 

Triangular 

Rectangular 

Rectangular 

Rectangular 

Triangular 

Rectangular 

Rectangular 


19.87 

0.0947 

Rectangular 

18.28 

0.3263 

Rectangular 

16.13 

0.6085 

Rectangular 

13.56 

0.8382 

Rectangular 

11.89 

0.8941 

Rectangular 

8.30 

1.2778 

Rectangular 

7.28 

1.6047 

Rectangular 

30.40 

23.20 


0.1170 
C . 1330 


Rectangular 

Rectangular 
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points  were  also  obtained  with  a  gas  gun,  but  the  target-to-f Iyer  thickness  ratio 
was  4:1  as  it  was  in  the  Boeing  experiments.  The  ETI  points  were  obtained  with 
an  exploding- foil  apparatus  under  the  PREDIX  program.  The  thickness  ration  was 
25:1  and  the  pulses  were  essentially  triangular.  McDonnell  Douglas  Astronautics 
Company  used  an  exploding  foil  apparatus  and  a  variable  thickness  ratio.  The 
last  two  points  were  obtained  by  Ferdman  and  Jajosky  also  with  an  exploding  foil. 
Although  some  of  these  results  are  relatively  old  and  different  definitions  of 
incipient  spall  were  used,  all  data  points  were  assigned  equal  weights.  No 
attempt  was  made  to  provide  more  recent  data  with  greater  significance  in  the 
analyses.  The  P-PUFF  computer  code  was  used  to  reduce  the  data  to  peak  tensile 
stress  and  pulse  duration.  Representative  stress-time  histories  for  each  data 
set  are  shown  in  figure  17  through  25.  Figure  26  is  a  plot  of  the  spall  thresh¬ 
old  showing  the  large  experimental  scatter  especially  for  very  short  pulse  dura¬ 
tions.  Each  of  the  models  was  again  fit  to  the  data  with  the  results  shown  in 
table  IX  and  figures  27  through  28. 

The  results  of  this  analysis  are  not .surprising.  Because  of  the  scatter  in 
the  available  data,  all  the  models  fit  the  data  reasonably  well.  This  is  obvious 
from  the  variances  listed  in  table  IX  and  from  an  examinatic  i  of  figure  27.  More 
discussion  of  this  point  is  presented  in  the  final  section. 

3.  FRACTURE  CALCULATIONS 

As  previously  mentioned  in  this  report,  the  goal  of  the  impact  experiments 
was  to  obtain  spall  data  on  6061-T6  aluminum  which  could  be  used  to  evaluate 
various  mathematical  models  of  dynamic  fracture  for  inclusion  in  PUFF-type  com¬ 
puter  codes  to  enhance  their  predictive  capabilities.  At  the  time  of  these 
experiments,  the  PUFF  fracture  criterion  was  critical  tensile  stress,  which  has 
been  shown  to  be  an  unrealistic  criterion  for  materials  of  interest. 

To  determine  the  adequacy  of  the  models  considered  for  predicting  dynamic 
fracture,  it  is  necessary  to  incorporate  them  into  a  PUFF-type  code  and  compare 
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Figure  17.  Stress-Time  History  at  Spall  Plane  (GM) 
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igure  Z*..  Scress-Time  History  at  Spall  Plane  (ETI) 


25.00 


L 

8 

in 

N 


© 

© 

m 


1 


S 


J- 

s 

» 

in 

i 


(AM)  3«n883«d 


o 

o 

in' 

T 


4! 


o 


-25.00 


30.00 


AFWL-TR-70-180 


TIME  [ftit cl 

25.  Stress-Time  History  at  Spall  Plane  (AFWL) 


AFWL  6 


! 


AFWL-TR- 70-180 


43 


At  (/tsec) 

Figure  26.  6061-T6  Aluminum  Spall  Data 
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Table  IX 


SPALL  MODEL 

CONSTANTS 

Model 

Assumption 

Variance 

Tuler-Butcher 

a.  (o  -  10. 206) t  =  1.069 

X  -  1.0 

13.0844 

b.  (o  -  5.487)t  -  4.979 

X  -  2.0 

12.1187 

c.  (o  -  4 . 250)  2 ' 2?t  *=  6.154 

o0  *»  4.25 

12.0973 

Rate  Process 

a.  a  =  10.011  -  10.906  logt 

12.6634 

Series  Expansion 

a.  o  -  9.255  +  1,557  _  0.034 

t 

12.5860 

experimental  results  with  results  calculated  using  each  of  the  models.  Thus, 
the  first  three  models  described  in  t*.ble  V  were  substituted  for  critical  ten¬ 
sile  stress  in  the  AFWL  P-PUFF  code.  The  changes  required  in  the  code  are  pre¬ 
sented  in  appendix  II.  A  version  of  SRI  PUFF  3  with  the  Hole-Growth  model  as  a 
fracture  criterion  was  available  and  was  used  for  the  evaluation  of  this  model. 

Sample  runs  were  made  for  each  model  using  0.020  inch  of  aluminum  impacting 
0.040  inch  of  aluminum.  Two  impact  velocities  were  used  for  each  model  such  that 
results  should  have  included  both  no-damage  and  incipient  spall.  According  to 
AFWL  experimental  results,  an  impact  velocity  of  0.20  mm/psec  should  not  have 
resulted  in  any  damage  while  a  velocity  of  0.26  mm/usec  was  required  for  incipient 
spall . 

Before  discussing  the  results  of  the  computations,  It  Bhould  be  noted  that 
some  difficulties  were  encountered  in  attempting  to  compare  calculations  based 
on  the  first  three  models  with  thoae  based  on  the  Hole-Growth  model.  These 
difficulties  arose  from  inherently  different  descriptions  of  fracture. 
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PUFF  and  P-PUFF  are  one-dimensional  Lagrangian  hydrodynamic  codes  which 
calculate  propagation  of  stress  waves  through  a  layered  homogeneous  material 
using  finite  difference  methods.  The  existing  fracture  routine  handles  spall 
by  creating  and  deleting  free  surfaces  at  zone  boundaries  when  certain  criteria 
are  met.  There  are  no  real  provisions  for  levels  of  damage  which  represent  less 
than  complete  separation  of  the  material.  Obviously,  this  is  not  physically 
realistic  since  damage  may  range  from  a  single  microcrack  to  complete  separation. 

The  Hole-Growth  model,  however,  is  capable  of  describing  various  levels  of 
damage  within  a  given  zone  by  either  increasing  void  density  or  increasing  void 
volume.  Thus,  in  discussing  the  results  of  fracture  computations,  it  is  appro¬ 
priate  to  postpone  those  for  the  SRI  Hole-Growth  model  until  later  In  the  section. 
The  other  three  models  are  similar  in  format  and  lend  themselves  to  simultaneous 
discussion. 

Table  X  presents  the  results  of  sample  runs  using  two  versions  of  the  Tuler- 
Butcher  model  as  well  as  the  Rate  Process  and  Series  Expansion  models.  All  four 
models  used  resulted  in  some  spall  even  at  the  low  velocity.  The  Tuler-Butcher 
model  appeared  to  be  relatively  Insensitive  to  impact  velocity  allowing  the  low¬ 
est  velocity  impact  to  spall  as  severely  as  the  incipient  shot.  The  Rate  Process 
and  Series  Expansion  models  were  better,  but  not  good. 

The  problem  which  arose  in  the  interpretation  of  these  results  wr.s  the 
definition  of  level  of  damage  in  the  computer  calculations.  It  was  not  clear 
whether  one  spalled  zone  should  be  defined  as  below  incipient  spall,  incipient 
spall,  or  complete  spall.  However,  no  immediate  decision  was  made  on  this 
question  because  of  the  discouraging  results  of  the  computed  shots. 

To  test  the  sensitivity  of  the  calculations  to  errors  in  the  experimentally 
determined  constants,  additional  computer  runs  were  made  with  the  constants  of 
each  model  individually  varied  by  5  percent  in  such  a  way  as  to  decrease  the 
number  ot  spalled  zones  which  resulted.  The  Tuler-Butcher  model,  as  expected, 


47 


AFWL-TR- 70-180 


Table  X 


NUMBER  OF 

SPALLED  ZONES  COMPUTED 

BY  FRACTURE  MODELS* 

Model 

.20  mm/psec 

.26  mm/usec 

.31  mm/psec 

Tuler 

Butcher,  X  =  1.81 

19 

20 

— 

Tuler 

Butcher,  X  =  2.00 

22 

22 

— 

Rate 

Process 

20 

26 

— 

Series  Expansion 

13 

17 

27 

*Targets  all  contained  200  zones. 

was  most  sensitive  to  the  value  selected  for  X.  A  5  percent  change  in  any  other 
constant  did  not  appreciably  affect  the  computed  results  for  any  of  the  models. 
Thus,  the  conclusion  was  that  none  of  these  models  could  provide  much  more  than 
a  crude  approximation  of  dynamic  failure  in  homogeneous  materials  of  interest. 

The  Hole-Growth  model,  however,  provided  more  encouraging  results.  Experi¬ 
mentally  determined  values  for  the  fracture  parameters  were  not  available  for  the 
6061-T6  aluminum  at  the  time  of  this  work.  The  values  used  were  approximations 
based  on  values  determined  for  1145  aluminum  by  SRI  (Ref.  4).  The  SRI  analysis 
is  based  on  nucleation  and  growth  of  voids  where  the  nucleation  rate  is 

N  (a,  t)  =  B  (a  -  c>N)  (4) 

and  the  growth  rate  is 

R  (o,  t)  =  C  (o  -  oj  R  (5) 

A  comparison  of  the  behavior  of  6061-T6  aluminum  undergoing  dynamic  failure 
with  that  of  1145  aluminum  reveals  a  number  of  characteristics  useful  in  approxi¬ 
mating  the  unknown  fracture  parameters.  Both  materials  undergo  ductile  dynamic 
failure  indicating  that  the  same  failure  models  are  applicable.  Studies  of  1145 
aluminum  indicated  that  voids  were  nucleated  at  inclusion  sites,  but  that  the 
number  of  potential  sites  was  far  greater  than  the  number  of  voids  observed.  Thus, 
even  though  the  6061-T6  aluminum  provides  a  still  greater  number  of  potential 
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nucleation  sites,  no  appreciable  change  in  the  value  of  B  in  equation  (4)  should 
be  required. 

However,  the  threshold  stress  for  nucleation,  ,  should  be  somewhat  higher 
than  that  for  1145  aluminum.  6061-T6  aluminum  does  require  slightly  higher 
stresses  to  cause  a  given  level  of  damage.  This  is  due  to  its  greater  yield 
strength  which  will  affect  the  threshold  stress  for  both  nucleation  and  growth. 

Finally,  growth  is  related  to  both  yield  strength  and  the  viscosity  of  a 
material.  The  growth  coefficient  is  inversely  proportional  to  the  viscosity  r, . 
Mitchell,  Hauser  and  Dorn  (Ref.  24)  have  shown  that  n  is  higher  for  6061-T6  alumi¬ 
num.  Thus  the  growth  coefficient  C  should  be  somewhat  lower. 

The  values  used  in  the  computer  calculations  are  summarized  in  table  XI.  The 
results  of  these  computations  are  illustrated  in  figure  29.  This  figure  also 
shows  photomicrographs  of  0.040-inch  targets  impacted  at  the  velocities  used  in 
the  calculations. 

Table  XI 

FRACTURE  PARAMETERS  FOR  6061-T6  ALUMINUM 

Zi 

B  —  5.0  x  10  No . /dyne-cm-sec 

S  2 

-  8.0  x  10  dynes /cm 

-3  2 

C  -  1.5  x  10  cm  /dyne-sec 

9  2 

C'q  -  4.0  x  10  dynes /cm 

From  figure  29,  it  is  apparent  that  the  threshold  for  incipient  spall  should 
be  defined  as  a  void  volume  of  approximately  2.5  percent  of  original  material 
volume.  To  further  verify  this  observation,  additional  calculations  were  per¬ 
formed  for  each  of  the  AFWL  data  points.  Results  are  shown  in  figure  30  with 
photomicrographs  of  actual  damage.  Since  all  damage  levels  are  approximately 
incipient  spall,  the  threshold  for  this  level  of  damage  might  be  defined  as  void 
volumes  of  2.2  to  3.5  percent  for  this  material. 
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0.20  mm/sec 


Figure  29.  Damage  Levels  Calculated  with 
Hole-Growth  Model  for  40  Mil  Target 
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But  perhaps  the  most  important  observation  to  be  drawn  from  these  comparisons 
is  the  ease  with  which  the  output  from  the  Hole-Growth  model  can  be  interpreted. 
The  other  models  considered  predict  inaccurate  and  physically  unrealistic  levels 
of  damage.  The  Hole-Growth  model  provides  a  numerical  value  for  the  level  of 
damage  based  on  the  relative  volume  of  voids  which  can  be  easily  correlated  with 
photomicrographs  of  actual  damage.  Damage  predictions  then  will  provide  a  real¬ 
istic  picture  of  the  actual  material  behavior. 
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SECTION  V 
CONCLUSIONS 

This  section  summarizes  the  results  obtained,  presents  the  most  obvious 
conclusions,  and  suggests  directions  for  future  research. 

1.  SUMMARY 

Three  separate  series  of  impact  experiments  were  performed  on  6061-T6  alumi¬ 
num  targets  on  the  AFWL  2.5  in.  gas  gun.  The  purpose  of  these  experiments  was 
to  determine  the  spall  threshold  for  the  material.  This  data  would  be  useful  in 
validating  various  models  of  dynamic  fracture  for  use  in  the  PUFF  computer  code. 

The  first  series  of  shots  showed  the  time-dependency  of  the  spall  threshold 
but  failed  to  bracket  the  threshold  for  the  thin  targets.  Certain  problems  with 
the  experimental  configuration  were  indicated,  and  a  second  series  of  shots  was 
planned . 

The  second  series  successfully  bracketed  the  spall  threshold.  A  comparison 
of  AFWL  data  with  that  generated  at  other  laboratories,  however,  revealed  dis¬ 
crepancies,  especially  for  short  duration  pulses. 

To  resolve  these  discrepancies,  a  third  and  final  series  of  shots  was  planned. 
6061-T6  aluminum  used  by  General  Motors  to  determine  a  spall  threshold  was  obtained 
for  the  AFWL  shots.  All  experimental  conditions  under  which  the  AFWL  data  was 
generated  were  maintained  so  that  the  material  was  the  only  variable.  The  results 
of  this  series  confirmed  the  GM  experiments  indicating  that  batch-to-batch  varia¬ 
tions  in  a  "standard"  commercially  obtainable  material  may  be  sufficient  to  cause 
widely  varying  results  in  dynamic  behavior. 

Experimental  results  were  evaluated  with  the  use  of  four  models  of  dynamic 
fracture:  the  Tuler-Butcher  model,  the  Rate  Process  model,  a  model  based  on  an 

inverse  series  expansion,  and  the  Hole-Growth  model  developed  by  SRI.  For  data 
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lrom  a  single  source,  the  Series  Expansion  model  provided  the  best  fit.  For  a 
large  number  of  data  points  from  a  variety  of  sources,  experimental  scatter  was 
severe  enough  that  no  model  provided  an  extraordinary  fit  and  all  were  acceptable. 

With  the  models  incorporated  in  PUFF-type  computer  codes,  calculations  were 
performed  to  determine  how  well  the  models  "predicted"  damage.  For  all  but  the 
Hole-Growth  model,  the  results  were  discouraging.  In  addition,  interpretation 
of  the  results  was  extremely  difficult.  The  Hole-Growth  model,  however,  provided 
easily  understandable  results  and  acceptable  accuracy  considering  the  approxima¬ 
tions  required. 

2.  CONCLUSIONS 

There  are  several  conclusions  to  be  drawn  from  this  effort.  First,  in  the 
study  of  dynamic  fracture,  careful  documentation  of  material  characteristics  and 
properties  as  well  as  interpretation  of  the  level  of  damage  is  required,  Batch- 
eo-batch  variations  in  the  properties  of  "standard"  6061-T6  aluminum  ’ed  to  a 
difference  in  spall  threshold  at  short  pulse  durations  of  approximately  15  per¬ 
cent  between  that  determined  by  AFWL  and  by  General  Motors.  Also,  while  much 
spall  data  for  6061-T6  aluminum  already  exists,  comparison  of  this  data  is  made 
difficult  by  lack  of  complete  definition  of  so-called  "incipient"  spall. 

The  data  analysis  presented  in  this  report  indicates  that  the  Hole-Growth 
fracture  model  developed  by  SRI  is  greatly  superior  to  other  existing  fracture 
models,  both  because  of  its  accuracy  and  because  of  the  ease  of  interpretation 
of  calculations. 

3.  RECOMMENDATIONS 

Because  of  the  superiority  of  the  SRI  model,  some  effort  should  be  made  to 
acquire  data  sufficient  to  obtain  the  fracture  parameters  for  at  lease  a  minimum 
number  of  materials  of  interest.  Proper  selection  of  materials  would  provide 
guidelines  for  classes  of  similar  materials.  SRI  has  already  used  data  on  1145 
aluminum  to  predict  results  of  experiments  on  copper.  The  predictions  were  within 
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10  percent  of  the  results  obtained.  It  Is  also  recommended  that  some  effort  be 
exerted  to  extend  the  model  to  include  effects  of  elevated  temperatures. 

Finally,  SRI  is  presently  extending  the  model  to  brittle  failure  in  Armco 
iron  and  graphite.  Thus,  satisfactory  failure  criteria  for  homogeneous  materials 
will  be  available.  Some  attempt  should  be  made  then  to  use  the  model  for  compos¬ 
ite  materials  of  interest.  These  materials  represent  the  real  problem,  and  no 
acceptable  failure  models  exist. 
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Similarly,  one  may  assume  a  value  for  cQ.  Values  for  the  other  constants  may 
then  be  determined  through  an  iterative  technique  which  evaluate  the  coefficient 
of  the  it  term  for  a  value  of  A.  A  standard  deviation  is  then  determined  for  each 
pair  of  constants  in  the  iteration.  The  pair  which  produces  the  minimum  standard 
deviation  is  the  solution.  The  rate  process  model 

o-a  +A0log(At)  (12) 

utilizes  only  two  constants  which  can  be  determined  directly  from  equations  (7) 
and  (8) . 

The  Series  Expansion  Model 

a  -  aQ  +  A/At  +  B/ (At) 2  (13) 

requires  evaluation  of  three  constants.  A  somewhat  more  complex  analysis  was  used 
for  this  model.  The  least-squares  mdethod  was  used  to  determine  all  three  con¬ 
stants.  The  pertinent  equations  are 


(14) 


(15) 


DENOM  =  n(T2T,-T32)  -  T,  (t^-T^)  +  T^Tj-T^) 


So  (T2T4-T32) 

1  '  T: 

(SjT.-SjTj)*  I2| 

[S1T3-S2T2) 

n  (S1Tu~S2T3) 

-  s0( 

DENOM 

t1t4-t2t3)  +  t2( 

.S2T1-S1T2) 

n(s2T;-S,T3; 

)-  T:l 

DENOM 

kS2T1-siT2)  +  S0 

(T1T3-V) 

DENOM 

(16) 


(17) 


(18) 


(19) 
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APPENDIX  II 

INCLUSION  OF  FRACTURE  MODELS  IN  P-PUFF  COMPUTER  CODE 

Each  of  the  fracture  models  discussed  in  the  tert  was  included  in  the  P-PUFF 
computer  code  in  order  to  determine  how  well  the  individual  models  "predicted" 
fracture.  The  original  version  of  P-PUFF  checked  for  fracture  in  the  HYDRO  sub¬ 
routine  if  any  zone  experienced  tensile  (or  negative)  stress.  The  failure  crite¬ 
rion  was  a  critical  tensile  stress  which,  if  exceeded,  caused  the  zone  to  "spall". 

To  replace  the  criterion  with  one  based  on  the  models  discussed,  it  is  neces¬ 
sary  to  compute  for  each  zone  and  each  cycle  some  function  of  stress  and  time 
whose  magnitude  can  be  compared  with  a  failure  constant. 

The  Tuler-Butclier  model  for  a  rectangular  pulse  was  used  by  determining 

whether  or  not  the  tensile  stress  in  a  zone  exceeded  the  threshold  stress  a  . 

o 

This  is  illustrated  in  the  flowchart  in  figure  31.  If  the  magnitude  of  the  ten¬ 
sile  stress  is  great  enough  and  the  zone  has  not  already  spalled,  a  quantity 
expressing  cumulative  damage  is  computed  as  follows: 

K  =  K  +  ^o-c0^  *DTNH  (20) 

where  DTNH  is  the  time  increment  for  the  cycle.  This  cumulative  quantity  is  com¬ 
pared  with  a  constant  which  must  be  exceeded  for  failure  to  occur.  When  failure 
does  occur,  the  existing  routine  is  used  (Ref.  7). 

Similarly,  a  cumulative  quantity  based  on  the  Rate  Process  model  can  be  used 
as  a  failure  criterion.  In  this  case  the  quantity  is  simply  the  duration  of  the 
pulse  for  tensile  stress  greater  than  threshold  stress  for  each  zone.  The  failure 
criterion  is 

o  -  A*log  j  q  ( DELT ( J) ]  >  cQ  (21) 

where  DELT(J)  is  puise  duration  for  zone  J.  Again,  the  existing  failure  routine 
is  used  when  failure  occurs.  Obviously  though,  this  criterion  is  extremely 
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Is  tensile  stress 
in  zone  Y  greater 
than  threshold  stress 
for  material  M? 


Has  the  zone  already 
spalled? 


Calculate  cumulative 
damage . 


If  cumulative  damage 
is  great  enough  for 
failure,  go  to  failure 
routine . 


Figure  31.  Flow  Chart  for  P-PUFF  Fracture  Criteria 
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sensitive  to  pulse  shape  in  chat  variations  in  stress  amplitude  severely  affect 
the  calculated  results.  For  truly  rectangular  pulses,  the  consequences  of  this 
are  minimized.  However,  experimental  pulses  are  not  truly  rectangular.  Thus, 
as  a  failure  criteria,  this  model  must  be  used  with  extreme  caution. 

The  Series  Expansion  model  also  requires  calculation  of  a  quantity  expressing 
cumulative  damage.  In  this  case 

K  *  K  +  ^-o+cq) *DTNH2  +  (A*DTNH)  (22) 

represents  cumulative  damage  for  each  zone.  This  quantity  is  compared  with  B 
which  is  defined  in  equation  (3)  in  the  text  to  determine  whether  or  not  failure 
occurs  in  a  zone. 

The  calculations  based  on  the  SRI  Hole-Growth  model  were  made  with  the  SRIPUF4 
computer  code.  For  details  concerning  Che  fracture  criterion,  the  reader  is 
referred  to  AFWL-TR- 70-99 ,  Dynamic  Fracture  Criteria  of  Homogeneous  Materials 
which  provides  more  complete  discussion  than  is  possible  here. 
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